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^ 1 ) Label root voxel with integer 0; 

2) Construct a queue and line up the root in the queue; 

IM* 3) If( There is at least one voxel in the queue ) 

py Serve the voxel x on the top of the queue: 

Q For( each of x's 26-connected neighbor voxel y) { 

i If( y in the volume and has not been labeled yet ) { 

Line up the y in the queue; 



/* label the voxel y */ 
Set dist = 999999 ; 

Fori each of y's 26-connected neighbor voxel 2 ) { 

If( z in the volume and has been already labeled with an integer of n z ) { 

where d(y,z) is 10, 14 or 17 if the Euclidean distance between y and z is 
, respectively; 

} 

If( <#« >d t ) { 

label y with integer dist ; 
dist = d t ; 

} 

} 

} 

} 

x leaves the queue; 

} 

Else { 

end of calculating the distance map. 

} 



